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DEVELOPED TURBULENT FLOW OF SUPERCRITICAL WATER WITH 
VARIABLE FLUID PROPERTIES IN A SMOOTH TUBE 
By Robert G. Deissler and Maynard F. Taylor 


SUMMARY 

A previous analysis of turbulent flow and heat transfer for air 
with variable properties flowing in smooth tubes is generalized in order 
to make it applicable to supercritical water. The generalization is 
necessary because all the pertinent properties of supercritical water 
vary markedly with temperature. Calculated velocity and temperature 
distributions, as well as relations among Nusselt number, Reynolds num- 
ber, and friction factor, are presented. The effect of variation of 
fluid properties across the tube on the Nusselt number and friction 
factor correlations can be eliminated by evaluating the properties at 
a reference temperature which ,is a function of both the wall temperature 
and the ratio of wall-to-hulk temperatures . 


INTRODUCTION 

An extension and generalization of the analysis for air with vari- 
able properties given in reference 1 were made in order to make it appli- 
cable to super critical water. (The term "supercritical water" designates 
water at pressures above the critical pressure and temperatures near the 
critical temperature.) An important characteristic of supercritical 
water, insofar as it affects heat transfer and flow, is the very large 
variation with temperature of its physical properties (specific heat, 
density, viscosity, and thermal conductivity) . Plots of the variation 
of the property values of water with temperature at a pressure of 
5000 pounds per square inch is given in figures 1 and 2 for data from 
references 2 to 4. From these plots it is evident that the radial 
variation of properties for supercritical water flowing in a tube with 
heat transfer should have a considerable effect on the velocity and tem- 
perature distributions, with a resultant effect on the heat-transfer 
coefficients and friction factors . 
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Several other analyses of heat transfer to supercritical water, 
which used assumptions not made herein, have been made and are 
reported in reference 3. 


ANALYSIS 


When the velocities and temperatures in a smooth tube are to be 
obtained as functions of distance from the wall, the differential equa- 
tions for shear stress and heat transfer may be written as follows: 


T 




CD 


q. = 



dt 

- PS Vh d^ 


( 2 ) 


Equations (l) and ( 2 ) written in dimensionless form, become 

t _ du _p e_ du 

T 0 ^0 dy + p 0 <3y + 

p 0 


(3) 


_a_ _ ( jl _L_ + _£_ °p g fr \ at + 

'So W) -^O p 0 c p,0 w J dy + 

. \ po/ 

(.The symbols used in this report are defined in the appendix.) 


W 


Assumptions . - The following assumptions are made in the use of 
equations (1) and (2) for obtaining velocity and temperature distribu- 
tions with heat transfer: 


(1) The eddy diffusivities for momentum and heat transfer s and 

e h are equal. Previous analyses for flow in tubes based on this assump- 
tion yielded heat-transfer coefficients and friction factors that agree 
with experiment except at low Peclet numbers. At low Peclet numbers the 
ratio of eddy diffusivities appears to be a function of Peclet number 
(Pe = Re Pr), but for Peclet numbers above 10,000, which is the range of 
interest for this analysis, the diffusivities are nearly equal (ref. 5 ). 

( 2 ) The expressions for eddy diffusivity that are found to apply to 
flow of air with heat transfer with variable properties (ref. l) apply 
also to supercritical water. These expressions are 

e = n 2 uy 


(5) 
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for flow close to the wall, and the Karman relation 



for flow feir from the wall, where n and x are experimental constants 
having the values 0.109 and 0.36, respectively. In the expression for 
flow close to the wall, the effect of kinematic viscosity p/p on e 
has been neglected. A preliminary investigation (not reported herein) 
indicates that neglect of this effect of p/p causes a somewhat larger 
variation of Nusselt number with Prandtl number than is indicated exper- 
imentally. 

(3) The variations of the shear stress t and the heat transfer 
per unit area q across the tube have a negligible effect on the veloc- 
ity and temperature distributions (ref. 1, fig. 12) . 

Flow close to wall . - For obtaining the velocity and temperature 
distributions close to a smooth wall, the expression for e from equa- 
tion ( 5 ) is substituted into equations (3) and ( 4 ) to give, in terms of 
the dimensionless parameters u + and y + . 


h Av 


p ‘t? _2 + +\ dt' 


pr 0 P 0 C p,0 


n u y 


where assumptions (l), ( 2 ), and (3) have been used. 
By rearrangement, equations ( 7 ) and (8) become 



(9) 

( 10 ) 
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Equations (9) and (lO) can then he written in integral form as 


u = 



dy 


1L+ P- 
Po p 0 


2 + + 
n u y 


(n) 


and 



dy 


JL _i_ + E- 


2 + + 
n u y 


(12) 


ko Pr 0 p O ^,0 

For a given wall temperature, k/kg, p/p Q , p/p©, and c p/ c p o are 
f mictions of t/tQj and since t/tQ = 1 - fit + ^ p/Po\ 


<>,0^0*0, 


^ k/kQ, 


p/ Pq, p/pq, and Cp/c q nre functions of f3t + . For given values of 


J 0> 

Lpq CLU.il P , 

iteration, that is, assumed relations between 


r* p *0 

equations (ll) and (l2) can be solved simultaneously by 

u + and y + and between 
t + and y + axe substituted into the right-hand sides of the equations 
and new values of u + and t + are calculated by numerical integration. 
These new values are then substituted into the right-hand sides of the 
equations and the process is repeated until the values of u + and t + 
do not change appreciably. Equations (ll) and (l2) give the relations 
among u + , t + , and y + for various values of the heat-transfer param- 
eter p for flow close to the wall. 

Flow at a distance from wall . - By use of assumptions (l), (2), 
and ( 3 ) , equat ions (3 ) and (4) become , for flow at a distance from the 
wall. 


i4 


1 = 


i*0 dy + p 0 

(a 2 uVdy +2 ) 


r- 
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k /*0 , P c p 2 

(du + /dy + ) 3 

dt + 

Pr o Po‘ c p,o 

(d 2 u + /dy +2 ) 

<3y + 


(13) 
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In equations (l3) and (l4), the molecular-shear-stress and heat -transfer 
terms have been retained so that the slopes for the velocity and temper- 
ature distributions close to and far from the wall can be matched at 
y^*. From the analysis for air presented in reference 1, it was found 
that y-|_ + , which is the lowest value of y + for which the equations 

for flow at a distance from the wall apply, has a value of 16 when the 
molecular-shear-stress and heat-transfer terms are retained. Substitu- 
tion of v for du + /dy + in equation (13) and rearrangement of the 
terms give 


dv 

v 


-x 


du 


V I ~ (h/mq)y 

pIpq 


which written in integral form is 



Cl5) 


(16) 


Substitution of v for du + /dy + in equation (l4) and rearrangement of 
the terms give 



(17) 


Elimination of dv/v from equations (15) and (l7 ) , rearrangement, and 
integration result in 




du 


+ v 


\ Pr 0 °V,0 ^0 


;) 


( 18 ) 
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If relations among v, u + , and. t + are assumed, equations (16 ) and 
(18) can be solved simultaneously by iteration. After the relation 

among v, t , and u + is determined, the expression v = — — or 

dy + 



(19) 


can be integrated, and y + can be calculated. 

Nub selt number, Reynolds number, and friction factor . - By use of 
the definitions of the quantities involved, it can be shown (see ref. 6) 
that with the fluid properties evaluated at the wall tempera-cure, the 
Nubs el t number, Reynolds number, and friction factor are given by 


Nuq = 

2r 0 Pr 

(20) 

V 

Re 0 = 

2r oV 

(21) 

II 

o 

2 

u + 2 

(22) 


where, for variable c^ and p. 



+ 



and 
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The relations between u + , y 4 *, and t + are calculated from equations (ll), 
(IE), (16), (18), and (is} . 

The Nusselt number, Reynolds number, and friction factor with the 
fluid properties evaluated at the bulk temperature rather than at the 
wall temperature can be found by using the relation t^/tQ = 1 - pt^ 4 ! 

Nusselt numbers, Reynolds numbers, and friction factors with properties 
evaluated at an arbitrary temperature in the fluid t x can be obtained 

from t x /t 0 = x(l - tfc/to) + t^tQ. 


RESULTS AND DISCUSSION 
Physical Properties 

The physical properties used in this analysis are shown in figures 1 
and 2. All the properties used were taken from tables and equations 
given in reference 2, with the exception of the thermal conductivities 
and viscosities at the lower temperatures which were taken from refer- 
ence 4. Considerable uncertainty exists in the property values of super- 
critical water, especially for viscosity and thermal conductivity, so 
that the accuracy of the computed results may be affected to some extent 
by errors in those values. The discrepancies among the viscosity and 
conductivity data used by various investigators are shown by the differ- 
ence between the solid and dashed lines in figure 2. 


Predicted Temperature and Velocity Distribution 

The predicted velocity and temperature distributions for various 
values of p for a wbUL temperature of 1360° R and a pressure of 
5000 pounds per square inch are shown in figures 3 and 4, respectively. 
Since the molecular -shear-stress and heat -transfer terms were retained 
in equations (l3) and (l4) for flow at a distance from the wall, curves 
in figures 3 and 4 were plotted using y^ 4 " equal to 16 as was used in 
reference 1 for air when the molecular -shear-stress and heat-transfer 
terms were retained. The peculiar variation of the physical properties 
of supercritical water necessitates the calculation of velocity and tem- 
perature distributions for each wall temperature . In general, the vari- 
ation of the velocity distributions is more regular than that of the 
temperature distributions. The irregular variation of the temperature 
distribution is due to the unusual variation of specific heat with tem- 
perature. In equation (l8) the specific heat appears in the denominator. 



8 


NACA EM E53B17 


Nusselt Numbers 

In figures 5 and 6, the predicted Nusselt numbers for supercritical 
water are plotted against Reynolds numbers for various values of the 
heat -transfer parameter p for a pressure of 5000 pounds per square 
inch and a wall temperature of 1360° R. 

Figure 5 shows the Nusselt numbers plotted against the Reynolds 
numbers with all the physical properties evaluated at the wall tempera- 
ture . The Nusselt numbers and Reynolds numbers are calculated from 
equations (20) and (2l) . At a given Reynolds number, the Nusselt num- 
bers can be seen to increase in a complex manner as p increases. In 
figure 6, all the physical properties in the Reynolds numbers and 
Nusselt numbers are evaluated at the bulk temperature. For a given 
Reynolds number, the Nusselt numbers decrease Irregularly with an 
increase in p. Because of the irregular variation of the Nusselt num- 
ber with P, no single reference temperature t can be found for the 
evaluation of the physical properties which will eliminate the effects 
of p. The effects of p can be eliminated, however, by using a ref- 
erence temperature which is a function of both the wall temperature and 
the ratio of wall -to -bulk temperature. In figure 7 the Nusselt numbers 
are plotted against the Reynolds numbers with the physical properties, 
including density, evaluated at the reference temperature for a pressure 
of 5000 pounds per square inch and various wall temperatures. These are 
the curves for p = 0 for the various wall temperatures and can be used 
for P ^ 0 by using the reference temperatures from figure 8, where 
x is plotted against the ratio of wall-to-bulk temperature for the 
various wall temperatures . The quantity x, which defines the reference 
temperature t x = x(tQ - t^) + t^, is calculated from the curves in fig- 
ure 5 and the corresponding curves for other wall temperatures . A refer- 
ence temperature can be selected for each point on the curves for 
P ^ 0 which brings the point to the p = 0 curve. The quantity x 
varies widely with both wall temperature and ratio of wall-to-bulk 
temperature . 

The variation of the level of the curves for Nusselt numbers 
plotted against Reynolds number in figure 7 is caused by the variation 
in Prandtl number. These curves indicate that Nusselt number varies 

approximately as PrQ^ ' ^ for the range of Prandtl numbers Bhown, 
whereas experimental data for various fluids indicate a variation of 

Ptq 0 - 45 . 1 This difference is probably a result of neglecting the effect 


^"An examination of the data for Prandtl numbers in the range considered 
in this report, for example, in references 7 and 8, indicates that an 
exponent of 0.45 fits the data better than the commonly used exponent 
0.4, which is an average for a larger range of Prandtl numbers than con- 
sidered herein. 
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of kinematic viscosity on the eddy diffus ivity in the region very close 
to the -wall; this region becomes increasingly important as the Prandtl 
number is increased. If a more accurate answer is desired, the present 
recommended procedure is to calculate Nu^ by multiplying the Nu^ ' s 

for a Prandtl number of 1 in figure 7 by the correction factor 

instead of using the curves for various wall temperatures in figure 7. 

The reference temperature for evaluating the properties in the Nusselt 
and Reynolds numbers can be obtained from figure 8. 

As an application of the foregoing analysis, wall and bulk temper- 
atures for flow through a 0.25 -inch tube with uniform heat input were 
calculated and are shown in figure 9. The bulk temperatures were obtained 
by setting the total heat transferred between the entrance and a given 
section equal to the change in enthalpy of the fluid passing the section 
and determining the temperature change corresponding to the enthalpy 
change (constant pressure) from the tables in reference 2. The curves 
for wall and bulk temperatures are nearly parallel at the lower temper- 
atures. As the bulk temperature approaches the critical temperature, 
the curves separate because of the rapid variation of physical properties 
in that region. The changes in slope of the bulk-temperature curve are 
caused by the very large increase in specific heat with temperature in 
the critical region and the subsequent decrease at higher temperatures . 

The wall -temperature curve does not follow the bulk-temperature curve 
because the properties of the fluid change from those of a liquid to 
those of a gas in passing through the critical region, with a consequent 
decrease in heat-transfer coefficient. It should be mentioned that the 
increase in specific heat near the critical temperature tends to increase 
the heat-transfer coefficient because of the larger Prandtl number In 
that regionj however, that effect is more than offset by the reduction 
in the values of the other properties . 


Several analyses which utilize assumptions other than those made in 
the present analysis were made and are reported in reference 3. Therein 
Goldman assumed that the curve of a plot of u + against y + for 3=0 can 


be used for flow with heat transfer if u + is redefined as / — 

*y i , — v-/o V T c/ p 


and y + Is written as 


f 

Jo 


^oA 


iVp 


d y, where local values of the prop- 


erties are used. Also, Elrod utilized Prandtl T s mixing- length 
theory and a laminar layer at the wall. The results of these two 
analyses are in reasonable agreement, as shown in reference 3. A com- 
parison between these analyses and the present one is made in figure 10. 
Most of the difference in calculated heat flux is caused by the differ- 
ence in the viscosities and thermal conductivities used (fig. 2). The 
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curves In figure 10 calculated from the present analysis and the prop- 
erties from reference 3 are in reasonable agreement with the analyses 
from reference 3. This agreement indicates that the calculated heat 
transfer is insensitive to many of the assumptions made in the analysis. 


Friction Factor 

The predicted friction factors plotted against Reynolds numbers for 
various values of the heat-transfer parameter p for a pressure of 
5000 pounds per square inch and a wall temperature of 1360° R are shown 
in figures 11 and 12. 

In figure 11, all the physical properties in the friction factors 
and Reynolds numb ers are evaluated at the wall temperature . For a given 
Reynolds number, an increase in the friction factors occurs with an 
increase in p. The friction factors plotted against Reynolds numbers 
with the physical properties evaluated at the bulk temperature are shown 
in figure 12. At a given Reynolds number, the friction factors can be 
seen to decrease nonuniformly with an increase in p . Since the friction 
factors vary in a complex manner with p, no single reference temperature 
t x can be found for the evaluation of the physical properties which will 
eliminate the effects of p. To eliminate the effects of p, a reference 
temperature which is a function of both the wall temperature and the 
ratio of wall-to-bulk temperature must be used. With all the physical 
properties evaluated at this reference temperature, the friction factors 
are plotted against the Reynolds numbers in figure 13. A plot of x 
against the ratio of wall-to-bulk. temperature for the various wall tem- 
peratures, where x is defined by t x = x(tg - t^) + t^, is shown in 

figure 14. The values of x were computed from the curves in figure 12 
and the corresponding curves for other wall temperatures. There is a 
variation of x with wall temperature and ratio of wall-to-bulk temper- 
ature for the friction factors, but it is not so great as the variation 
for the Nusselt numbers. 


SUMMARY OF RESULTS 

The following results were obtained from the analytical investigation 
of turbulent flow and heat transfer of supercritical water in smooth 
tubes : 

1. Heat addition to the fluid caused a flattening of the velocity 
and temperature profiles . The variation of temperature distribution 
with heat addition was complex because of the unusual variation of spe- 
cific heat with temperature. 
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2. The effects of variable properties on the Nusselt number and 
friction factor correlations were eliminated by evaluating the prop- 
erties of the Nusselt number, Reynolds number, and friction factor at 
reference temperatures which are functions of both wall temperature and 
ratio of wall -to -bulk temperature . 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX - SXMBOIS 

The following symbols axe used in this report: 

specific heat of fluid at constant pressure , Btu/(lb)(°F) 

specific heat of fluid at constant pressure evaluated at t n , 
Btu/ (lb } (°F ) 

inside diameter of tube, ft 

base of natural logarithms 

acceleration due to gravity, 32.2 ft/sec^ 

heat-transfer coefficient, r r— , Btu/(sec)(sq ft)(°F) 

^0 " Td 

thermal conductivity of fluid, Btu/(sec)(sq ftJC^F/ft) 

thermal conductivity of fluid evaluated at tu, 

Btu/(sec)(sq ft)(°F/ft) 

thermal conductivity of fluid evaluated at t Q , 

Btu/(sec)(sq ft)(°F/ft) 

thermal conductivity of fluid evaluated at t , 

Btu/(sec)(sq ft)(°F/ft) 

constant 

static pressure, lb/sq ft abs 

rate of heat transfer toward tube center per unit area, 
Btu/(sec)(sq ft) 

rate of heat transfer at inside wall toward tube center per 
unit area, Btu/(sec)(sq ft) 


Tq inside tube radius, ft 

t absolute static temperature, °R 

tQ absolute wall temperature, °R 

t bulk or average static temperature of fluid at cross section of 

tube , °R 

t film temperature, x(t Q - t^) + t^, °R 
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u time -average velocity parallel to axis of tube, ft/sec 

u b bulk or average velocity at cross section of tube, ft/sec 

x number used in evaluating arbitrary temperature in tube, t x 

y distance from tube wall, ft 

6 coefficient of eddy diffusivity for momentum, sq ft /sec 

coefficient of eddy diffusivity for heat, sq ft/sec 
x Karman constant 

p absolute viscosity of fluid, (lb)(sec)/sq ft 

p^ absolute viscosity of fluid evaluated at t^, (lb) (sec )/sq ft 

Pq absolute viscosity of fluid evaluated at tQ, (lb) (sec)/ sq ft 

P x absolute viscosity of fluid evaluated at t x , (lb)(sec)/sq ft 

p mass density, (lb)(sec^)/ft^ 

p.^ b ulk or average density at cross section of tube, (lb) (sec^)/ft- 
Pq mass density of fluid at wall, (lb ) (sec^ ) / ft 4 
P x density of fluid evaluated at t x , (lb)(sec^)/ft^ 

T shear stress in fluid, Ib/sq ft 

t q shear stress in fluid at wall, Ib/sq ft 

Subscripts: 

fr on friction-pressure gradient 


D(dp/dx) fr 2 t q 

t 0^ “ 2~ = 2 
2p 0 p 0 u b 


Dimensionless groups : 


JLVX p K vwyw • ■ 

wA 

heat -transfer parameter, — — 


friction factor with density evaluated at 
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2T r 


Pe 

Pr 

Pr C 

Re 

Re r 


■‘o 

t + 

*1 + 

V 

u + 


friction factor with density evaluated at tp. 


friction factor with density evaluated at t x . 


Pb u b 
2^0 

P X V 


Ruq Nusselt number with thermal conductivity evaluated at 

Nujj Nusselt number with thermal conductivity evaluated at t^, 

NUjj. Nusselt number with thermal conductivity evaluated at t x . 


Peclet number 
Prandtl number 

Prandtl number with properties evaluated at wall, c^ q p.Qg/kQ 
Reynolds number 

Reynolds number with density and viscosity evaluated at t Q , 

PqV* 

MO 

Ke^ Reynolds number with density and viscosity evaluated at t^, 

PbV 3 

Mb 

Re x Reynolds number with density and viscosity evaluated at t v , 


PxV 5 


x' 


VV S 0 

tube -radius parameter, \Iq/&q r ° 


temperature parameter > 


value of t + at y^ + 


(t 0 -t) Cp^Q gT 0 1 - t/t Q 


hulk -temperature parameter 

u 




velocity parameter. 


-\Ao/po 


>ris <ris <ri6 
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u^ bulk velocity parameter. 


u-L + value of u + at y^ 4 " 


^b 


7 ^ 7 ? 


0 


du 

dy 4 


value of v at 


*1 


"vAo/Pq 

wall- distance parameter, ^ 


value of y at intersection of curves for flow close to wall 
and at a distance from wall 
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Figure 1. - Variation of density and specific heat of water with temperature. 
Pressure, 5000 pounds per square inch. Data obtained from references 2 
and 4. 
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Figure 2. - Variation of thermal conductivity and viscosity of water with tempera- 
ture. Pressure, 5000 pounds per equate Inch. 
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Figure 3. - Predicted generelited Telocity distributions for flow of super critical 
5000 pounds per square Inch; vail temperature, 1360° R. 





* 



Figure 4. - Predicted generalized taupe ret ore distributions far flav of supercritical vater vlth heat addition. Pressure 
5000 pounds par square Inch; vail temperature, 1300° B. 










Figure 6. - Predicted variation of Nusselt number with Reynolds number for flow of supercritical 
water with heat addition and properties evaluated at bulk temperature. Pressure, 5000 pounds 
per square Inch; wall temperature, 1360° R. 











Distance along tube, in. 


Figure 9. - Predicted variation of vail and huUc temperatures for flow through 0, 25-inch tub© 
with constant heat input using curves in figure 7. Heat input, 575 Btu per second per 
square foot; weight flow, 0.0198 pound per second; pressure, 5000 pounds per square inch. 
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Figure 10. - Comparison of heat flux predicted by reference 3 and in present 
analysis for wall temperatures of 1260° and 1460° R. Tube diameter, 
0.00521 foot; weight flow, 0.0118 pound per second; pressure, 5000 pounds 
per square inch. 
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